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GALAXY EVOLUTION 


Isolated compact elliptical galaxies: 
Stellar systems that ran away 


Igor Chilingarian^’^* *and Ivan Zolotukhin® 


Compact elliptical galaxies form a rare class of stellar system (~30 presently known) 
characterized by high stellar densities and small sizes and often harboring metal-rich 
stars. They were thought to form through tidal stripping of massive progenitors until two 
isolated objects were discovered where massive galaxies performing the stripping could 
not be identified. By mining astronomical survey data, we have now found 195 compact 
elliptical galaxies in all types of environment. They all share similar dynamical and stellar 
population properties. Dynamical analysis for non-isolated galaxies demonstrates the 
feasibility of their ejection from host clusters and groups by three-body encounters, 
which is in agreement with numerical simulations. Hence, isolated compact elliptical and 
quiescent dwarf galaxies are tidally stripped systems that ran away from their hosts. 


G alaxies are thought to form through 
the hierarchical merging of smaller 
building blocks into larger systems 
{1, 2) and the history of these in¬ 
teractions is imprinted in their ob¬ 
servable properties. Some galaxies such as 
ultra-compact dwarfs {3) and compact ellipti¬ 
cals (cEs) {4-7), show evidence of strong tidal 
interactions with massive neighboring galax¬ 
ies (8) that stripped most of the stars from the 
compact galaxies’ progenitors. cEs are rare 
galaxies with high stellar densities that resem¬ 
ble centers of giant ellipticals but have masses 
that are about two orders of magnitude smaller 
[M ~ 10® solar mass (M©)]. They are found 
mostly in the cores of galaxy clusters next to 
massive central galaxies, which is in alignment 
with the above hypothesis for their evolution. 

The recent discoveries of isolated cE galax¬ 
ies (9, 70) that do not belong to any galaxy 
cluster or group raised another round of debate 
about cE formation: whether they all formed 
through the tidal stripping or through a differ¬ 
ent mechanism of formation such as mergers 
of dwarf galaxies with specific morphologies 
and configurations (70). Dwarf-dwarf galaxy 
mergers do happen in vicinities of massive 
galaixes (77, 72). However, neither have they 
been observed in low-density environments, 
nor any of the remnants resemble properties 
of cE galaxies. The existence of a substan¬ 
tial number of isolated cEs will hence imply 
notably higher dwarf-dwarf merger rates than 
predicted by numerical simulations (72) and 
challenge the currently accepted hierarchical 
structure formation paradigm. 

We demonstrated that all known cE galax¬ 
ies are outliers from the universal optical- 
ultraviolet color-color-magnitude relation of 
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galaxies (74). We could therefore perform a 
search for cE galaxies not only in the centers 
of rich clusters and groups as has been done 
before (5), but across all environments using 
data from wide-field imaging surveys, the op¬ 
tical ground based Sloan Digital Sky Survey 
[SDSS, (75)] and the ultra-violet all-sky survey 
carried out by the GALaxy Evolution eXplorer 
[GALEX, (76)] spacecraft, which are all pub¬ 
licly available in the Virtual Observatory. 

First, we created an initial list of candi¬ 
dates (supplementary materials) from the sam¬ 
ple of galaxies having spectra in the SDSS 
and, hence, known distances, by selecting out¬ 
liers above -i-0.035 mag in the optical (g — r) 
color from the universal relation (74). We 
chose low luminosity galaxies [L < 4 x 10® 
solar luminosity (7)©) or absolute magnitude 
(Mg) > —18.7 mag] that had small half-light 
radii {Re < 0.6 kpc) or were spatially unre¬ 
solved in SDSS images; did not show substan¬ 
tial ellipticity, which was essential for remov¬ 
ing edge-on spiral galaxies; had the redshifts in 
the range of 0.007 < z < 0.08 (distances be¬ 
tween 30 and 340 Mpc); and either possessed 
red near-ultraviolet colors [{NUV — r) > 
4 mag] or remained undetected in the NUV 
band. We constrained by color and also re¬ 
moved objects that have emission lines in their 
spectra in order to exclude any objects with re¬ 
cent or ongoing star formation. 

We then fitted their SDSS spectra against 
a grid of stellar population models using the 
“NBursts” code (77) and obtained mean ages, 
metallicities, and velocity dispersions of their 
stars. We rejected candidates with stellar ages 
younger than 4 billion years and introduced an 
additional constraint based on stellar velocity 
dispersions (a > 60 km/s). Stellar systems 
in equilibrium that are dynamically supported 
by random motions of stars, as most ellipti¬ 
cal galaxies are have their dynamical masses 
(Mvir), half-mass radii {Re), and global ve¬ 
locity dispersions ((j„) connected by the sim¬ 


ple virial relation: Mvir = 9.96i?e(r®/G 
{18, 19). Therefore, for a galaxy with known 
velocity dispersion and a stellar mass (M*) de¬ 
rived from its luminosity and stellar population 
parameters, we can estimate the lower limit 
for the half-light radius (if a galaxy contains 
dark matter, its real half-light radius will be 
larger because Mvir > M*). Hence, we can 
firmly reject physically extended objects such 
as “normal” dwarf elliptical galaxies that are 
unresolved in SDSS images because of their 
large distances by selecting only objects with 
high velocity dispersions. 

In this fashion, we constructed a sam¬ 
ple of 195 galaxies (Fig. 1 and supplemen¬ 
tary materials). We then cross-matched this 
list against the SDSS Galaxy Groups catalog 
(20) and established their group/cluster mem¬ 
bership. For seven objects without counter¬ 
parts in the Galaxy Groups catalog, we iden¬ 
tified possible host galaxies (in most cases, 
group centers) located between 750 kpc and 
3.3 Mpc in projection. Because some bright 
and extended galaxies were missing from the 
SDSS spectroscopic sample, and therefore also 
from the Galaxy Groups catalogue, we used 
NASvVInfrared Processing and Analysis Cen¬ 
ter (IPAC) Extragalactic Database (NED) for 
the identification of host galaxies for 45 cEs. 
Our final sample contains 56 galaxies in clus¬ 
ters, 128 in groups, and 11 isolated or field cE 
galaxies. Eight galaxies (supplementary mate¬ 
rials) exhibit prominent tidal streams similar to 
the two known cEs (8). 

Ages and metallicities (Fig. 2 and sup¬ 
plementary materials), colors and luminosi¬ 
ties of 11 isolated cE galaxies do not show a 
statistically significant difference from those 
of galaxies being members of clusters and 
groups. The Kolmogorov-Smimov (KS) prob¬ 
abilities of these properties for isolated and 
non-isolated subsamples to be derived from 
the same underlying distribution range from 30 
(for luminosities) to 80% (metallicities). Our 
cE sample does not show any correlation be¬ 
tween the metallicity and the stellar mass con¬ 
versely to normal elliptical galaxies, which ex¬ 
hibit a rather tight mass-metallicity relation 
(Fig. 2). 

For rich groups and clusters hosting cE 
galaxies, we built caustic diagrams (21-23) 
that present differences of radial velocities of 
cluster members from the cluster center versus 
projected distances. A galaxy position on such 
a diagram reflects its dynamical status: objects 
sitting deep inside the cluster potential well are 
located inside the distribution whereas galaxies 
near the edges are barely gravitationally bound 
to the host cluster or infalling on to it for the 
first time. 

We constructed an ensemble cluster by nor¬ 
malizing individual cluster and group data by 
their velocity dispersions and sizes for 33 
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Fig. 1: A sample of 
compact elliptical 
galaxies in all types of 
environment. Our current 
sample of compact 
elliptical galaxies (blue 
symbols) is compared with 
a dataset compiled from 
the literature (red symbols). 
Dots and stars denote 
group/cluster and isolated 
compact elliptical galaxies 
respectively. Square 
panels in the bottom part of 
the figure show 
representatives of the 
current sample, and top 
corner insets display 
Messier 32 and NGC 4486, 
prototypical cEs in the local 
Universe as they would 
look with the SDSS 
telescope at a 130 Mpc 
distance {z = 0.03). Each 
inset panel covers a 20 by 
20 kpc region centered on 
a cE. 
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Fig. 2: Comparison of 
stellar population 
properties of isolated 
and nonisolated 
compact elliptical 
galaxies to a 
reference sample of 
elliptical galaxies 
from the SDSS. Mean 
stellar ages and 
metallicities of cE 
galaxies hosted by 
groups and clusters 
(round symbols) and 
isolated cEs (stars) 
were determined from 
the modeling of their 
SDSS spectra by using 
the “NBursts” full 
spectrum fitting. The 
plotting symbols 

are color coded according to the r-band absolute magnitudes derived from the SDSS photometry. 
We also computed ages and metallicities for a reference sample of 82500 elliptical galaxies from 
SDSS DR7 in the same fashion along with the median r-band absolute magnitude for every bin of 
age-metallicity parameter space. These magnitudes are shown as a background map, with the 
colors on the same scale as cE galaxies. This figure demonstrates that: (i) stellar populations of 
isolated and group/cluster cEs do not differ statistically and (ii) cEs are on average much fainter 
than are normal elliptical galaxies of the same stellar age/metallicity. 
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structures from our sample each of which in¬ 
cluded over 20 member galaxies (24, 25). Then 
we computed its caustic diagram in order to 
visualize the phase space pattern of the in¬ 
falling galaxy population and overplotted our 
cE galaxies on it (Fig. 3). The cE popula¬ 
tion strongly differs from other cluster mem¬ 
bers. The KS tests for projected distance and 
radial velocity distributions reject the hypoth¬ 
esis of cE and cluster member samples being 
derived from the same parent population at the 
97 and 98% levels. Numerical simulations of 
tidal stripping (5, 26) suggest that a progenitor 
galaxy, even if it approaches a cluster center on 
a very extended radial orbit, will lose a major 
fraction of its orbital energy because of dynam¬ 
ical friction, become gravitationally locked in 
the inner region of a cluster on a tightly bound 
orbit, and will finally be accreted by the host 
galaxy after a few billion years. Many cEs 
from our sample exhibit this behaviour (Fig. 3 
and supplementary materials). However, we 
see a number of cE galaxies close to the edges 
of the caustics suggesting that they are barely 
gravitationally bound to the cluster potential 
because they do not belong to the infalling pop¬ 
ulation as we demonstrated. This looks com¬ 
pletely unrealistic in the case of a one-to-one 
galaxy encounter resulting in tidal stripping, 
but in the case of a three- or multiple-body 
encounter this situation becomes significantly 
more likely. 

An interaction of binary stars with the cen¬ 


tral supermassive black hole is one accepted 
scenario for the creation of hypervelocity stars 
(27) in our Galaxy: One of the binary com¬ 
ponents is ejected, whereas the other one falls 
on-to the black hole. Numerical simulations 


suggest (28) that three-body encounters are re¬ 
sponsible for putting Milky Way satellites on 
extreme orbits going as far as 3 Mpc away. 
Even though typical galaxy clusters have much 
wider and deeper potential wells than that of 
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Fig. 3: Positions of 
cE gaiaxies identified 
on caustic diagrams 
in gaiaxy clusters 
and rich groups with 
more than 20 
members in the 
SDSS. (Top) An 
example of a caustic 
diagram showing 
projected distances 
and radial velocity 
differences of the 
members of the cluster 
MKW 03s that includes 
cE galaxies shown as 
stars. The derived 
caustics shown as 
solid blue lines roughly 
correspond to the 
escape velocity for a 
galaxy at a given 
distance from the 
cluster center. Shaded 
areas show statistical 
1 cr uncertainties of the 
caustic line 

computation. (Bottom) 

The caustic lines for an 
ensemble cluster of 
2,592 individual 
members are shown as 
gray shaded areas and 
light gray dots, 
respectively. The 
caustic lines 
normalized by 

corresponding velocity dispersions and half-mass radii (Rgff) are shown for six individual clusters 
as colored solid lines. Small stars denote cEs located deeply inside the potential wells of these 
six clusters, and large stars (with same colors as caustic lines) indicate cEs that have projected 
radial velocities of at least 85% of the caustic amplitude. Vectors point to the location of identified 
host galaxies for corresponding cEs. Three of them are different than the cluster central galaxies 
indicating that these cEs belong to subgroups inside clusters and that the caustic diagnostic is 
irrelevant for them. However, the four circled cEs are examples of galaxies barely gravitationally 
bound to their host clusters. 
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the Local Group, three- and even multiple- 
body encounters must happen much more fre¬ 
quently in those dense environments. There¬ 
fore a certain probability exists that some of 
them will lead to the gravitational ejection of 
galaxies participating in the interaction to ex¬ 
treme radial orbits with the apocentric dis¬ 
tances of a few megaparsecs (29). A three- 
body encounter that might eject a cE galaxy 
from its host cluster or group does not have to 
happen during the cE formation through tidal 
stripping, that is two galaxies do not have to 
fall on-to the cluster/group center at the same 
time. When a cE progenitor is tidally stripped, 
it quickly settles on a tightly-bound rapidly 
decaying orbit (5) and if another galaxy in¬ 
falls later, but before a newly formed cE has 
been accreted (hundreds of millions to a cou¬ 
ple of billion years), the three-body encounter 
becomes possible. 

We estimate the probability of a close three- 
body encounter geometrically. Numerical sim¬ 
ulations suggest {30) that over a typical cE life¬ 
time of 2 billion years (5), an average brightest 
cluster galaxy (BCG) must have experienced 
three of four mergers with massive (M > 
10^° Mq) galaxies. We assume that: (i) a typ¬ 
ical cE resides on a quasicircular orbit within 
r^E ~120 kpc from a host BCG galaxy af¬ 
ter correction for projection effects (fig. S4), 
(ii) galaxies infall on a BCG on radial or¬ 
bits from random directions, and (iii) a three- 
body encounter will be sufficiently close if a 
cE passes within rat « 20 kpc from a mas¬ 
sive infalling galaxy. Hence, the probability 
is as a volume ratio of a cylinder of radius 
rat, height rcE and a sphere of radius r^E, 
P-ib = 3/4(rat/rcB)^ « 0.02, or ~6 to 8% 
for three of four merger events. 

In our sample of cluster and group cE 
galaxies, we indeed see numerous examples in 
which a cE resides only 20 to 80 kpc in pro¬ 
jection from an ongoing major merger scene 
or several other massive cluster members are 
visible in the cE vicinity apart from the mas¬ 
sive central cluster/group galaxy. Also, there 
is a known example of a globular cluster in the 
Virgo cluster {31) that was likely ejected at the 
speed of 2500 km/s and became gravitationally 
unbound to the cluster and its central galaxy 
Messier 87. 

We conclude that the tidal stripping process 
can explain all observational manifestations 
of compact elliptical galaxies, including the 
formation of isolated cEs whose existence was 
suggested as a strong counter-argument for 
tidal stripping (9). The ejection of cEs from 
central regions of galaxy clusters by three- 
body encounters is a channel for these galaxies 
to survive for an extended period of time in 
the violent cluster environment where they 
would otherwise be accreted by massive hosts 
on a timescale of 2 billion to 3 billion years. 


The 11 isolated cEs probably represent a pop¬ 
ulation of runaway galaxies that received suffi¬ 
cient kick velocities to leave their host clusters 
or groups forever. 

The gravitational ejection mechanism may 
also explain the very existance of extremely 

rare isolated quiescent dwarf galaxies {32), 
where the star formation quenching is usually 
explained by environmental effects. These sys¬ 
tems are more spatially extended than cEs and 
do not exhibit substantial tidal stripping foot¬ 
prints. This suggests that they never came very 
close to cluster/group centers and therefore, 
the three-body encounter probability for them 
should be lower than that for cEs although still 
nonnegligible. 
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SUPPLEMENTARY MATERIALS 

51 Input galaxy catalogue 

As a starting point for the search of cE galaxy 
candidates, we compile the input galaxy cat¬ 
alogue. In order to do this, we use pub¬ 
licly available datasets, GALEX Data Release 
6 (16), SDSS Data Release 7 (15) and UKIDSS 
Data Release 8 (33). Initially, from the SDSS 
DRV spectral sample we select normal non¬ 
active galaxies with the redshift in the range 
0.007 < z < 0.6. This list is cross-matched 
with the UKIDSS Large Area Survey cata¬ 
logue using the best match criterion with the 
3 arcsec radius. The SDSS-UKIDSS sample 
is then uploaded to the GALEX CasJobs ser¬ 
vice in order to look for best matches within 
3 arcsec. As a result of this procedure, we 
obtain a catalogue of 429 707 galaxies with 
fully corrected UV-to-optical spectral energy 
distributions in 11 photometrical bands (FUV 
and NUV from GALEX, ugriz from SDSS, 
and Y JHK from UKIDSS) and spectra. Pet¬ 
rosian magnitudes are fc-corrected using an¬ 
alytical approximations (34). SDSS spectra 
are fitted with state-of-the-art stellar popula¬ 
tion models. This procedure yields ages and 
metallicities of their stars as well as stellar ve¬ 
locity dispersions which we add to the cata¬ 
logue. 

52 Full spectral fitting 

We use the “nbursts” full spectral fitting tech¬ 
nique (17) with high-resolution (R = 10,000) 
PEGASE.hr (35) simple stellar population 
(SSP) models computed using the Salpeter 
stellar initial mass function (36) in order to ex¬ 
tract kinematics and stellar population proper¬ 
ties from our spectra (see Pig. [S3. 

S3.1 Initial sample 


The fitting procedure comprises the follow¬ 
ing steps: (1) a grid of SSP spectra with a 
fixed set of ages (spaced logarithmically from 
20 Myr to 18 Gyr) and metallicities (from —2.0 
to +0.5 dex) is convolved with the wavelength- 
dependent instrumental response of SDSS pro¬ 
vided along with every spectrum in the 6th 
FITS extension; (2) a non-linear least square fit¬ 
ting against an observed spectrum is done for 
a template picked from the pre-convolved SSP 
grid interpolating on age (log t), and metallic- 
ity (Z), broadened with the line-of-sight ve¬ 
locity distribution (LOSVD) parametrized by 
V and a and multiplied pixel-by-pixel by an 
order Legendre polynomial (multiplicative 
continuum), resulting in n + 5 parameters to be 
determined (we used n = 19 for SDSS spec¬ 
tra). 

The fitting procedure, the error analysis and 
degeneracies between kinematics and stellar 
population parameters have been throughly de¬ 
scribed (37). We analysed the systematics of 
stellar population parameters with respect to 
the a-element enhancement (38), the sensi¬ 
tivity of different spectral features (39). The 
choice of the stellar initial mass function in 
the models is shown not to affect stellar metal- 
licity estimates and only slightly affect stellar 
ages (40). 

In order to precisely estimate velocity dis¬ 
persion, we convolve the SSP grid with the 
spectral line spread function (LSE) of the spec¬ 
trograph varying across the wavelength range. 
In SDSS DR9 spectra, the LSE shape is avail¬ 
able in the 6th EITS extension of every data 
file. It was demonstrated (41) that for inter¬ 
mediate signal-to-noise spectra the stellar ve¬ 
locity dispersion can be estimated with 10% 
uncertainties down to ~l/2 of the LSE width 


(expressed as Gaussian a). For SDSS this cor¬ 
responds to 45 km/s that is substantially be¬ 
low our velocity dispersion selection criterion 
(a > 60 km/s, see below). 

S3 Selection of candidate compact elliptical 
galaxies 

By placing known cE galaxies on the 3D 
optical-UV color-color-magnitude diagram 
we noticed that they occupy a distinct region in 
this parameter space (14). We, hence, use sev¬ 
eral approaches to define this region in order 
to make it suitable for a formal database query. 
They result in 5 samples of cE candidates ex¬ 
tracted from the catalogue described above. 
Since the input data were generated by several 
sometimes imperfect data reduction pipelines, 
we have to visually inspect every candidate im¬ 
age in order to discard extended galaxies with 
bright nuclei, edge-on spiral galaxies and arte¬ 
facts, which also satisfied our initial automated 
cE selection criteria. 

We selected objects that have: (a) red opti¬ 
cal (ff — r) colors, at least +0.035 mag above 
the universal color-color-magnitude relation 
(14); (b) low luminosities [L < 4 • 10®I/Q or 
Mg > —18.7 mag]; (c) small half-light radii 
[i?e <0.6 kpc] or remain spatially unresolved 
in SDSS images; (d) no significant ellipticity; 

(e) the redshifts in the range 0.007 < a < 0.08 
[distances between 30 and 340 Megaparsecs]; 

(f) either possess red near-ultraviolet colors 
[{NUV — r) > 4 mag] or drop out of the 
NUV band; (g) no significant emission lines 
in their spectra; (h) velocity dispersions a > 
60 km/s and stellar population ages / > 4 Gyr. 

In the next subsections we list the criteria 
that we used to build those samples. 


The initial sample included 108 galaxies that have NUV — r >4.0 mag, the absolute g magnitude fainter than -18.7 mag, the 50 per cent Petrosian 
radius of smaller than 2 arcsec in the angular measure or below 0.7 kpc in the physical scale, a distance from the 3D surface (14) along the g — r 
color axis greater than 0.035 mag, the stellar age older than 4 Gyr and the velocity dispersion above 60.0 km s“^, the latter two parameters derived 
from the SSP fitting. 

Below we give the formal specification of this filter using subsets expressions from the TOPCAT table processing software (42) with the column 
names from the catalogue. Note that GR_FIT column is a 3D surface value in (g — r, NUV — r. Mg) coordinates. 

(corrmag_NUV - corrmag_r - kcorr_NUV + k;corr_r) > 4,0 && 

(corrmag_g - kcorr_g - 25 - 5 * loglO(luminosityDistance(z, 72.0, 0.3, 0.7))) > -18.7 && 

(petror50_r < 2.0 I I petror50_r / 206.265 * luminosityDistance( z , 72.0, 0.3, 0.7) < 0.7) && 

(corrmag_g - corrmag_r - kcorr_g + kcorr_r - gr_fit) > 0.03 && 
age_ssp > 4000.0 && dispvel_ssp > 60.0 

After visual inspection, we have rejected 52 contaminated objects and therefore 56 cE galaxies were left in this sample. The latter two criteria 
on the age and velocity dispersion (age_ssp > 4000.0 && dispvel.ssp > 60.0) are common between all samples except the sample 
selected by the virial radius criterion. 


S3.2 UV-dropout sample 

The second sample is constructed with a relaxed requirement of a candidate detection in the NUV band, because only 47 per cent of galaxies from 
our catalogue have NUV magnitudes. We dropped the constraints on the NUV — r color and the one derived from the 3D photometric relation 
(because it also requires a NUV magnitude). Instead, together with our usual requirement of a galaxy having SSP age older than 4 Gyr and the 
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velocity dispersion higher than 60 km s”*^, we added a simple geometrical constraint from the canonical color-magnitude {g — r, Mr) diagram: a 
candidate galaxy is redder than the red sequence plus la scatter. This translates into the following expression: g — r > 0.8 — (Mr + 20.0) ■ 0.1/3.0. 
Hence, we obtain the following formal query expression: 

(corrmag_g - kcorr_g - 25 - 5 * logl0(luminosityDistance(z, 12,Qf 0.3, 0.7))) > -18.7 && 

(petror50_r < 2.0 I I petror50_r / 206.265 * luminosityDistance{z, 72.0, 0.3, 0.7) < 0.7) && 
age_ssp > 4000.0 && (dispvel_ssp >60.0 I I dispvel_exp > 60.0) && 

(corrmag_NUV < -100 I I NULL_corrmag_NUV) && 

(corrmag_g - corrmag_r - kcorr_g + kcorr_r > 0.8 - 

(corrmag_r - kcorr_r - 25 - 5 * logl0(luminosityDistance(z, 72.0, 0.3, 0.7)) + 20.0) * 0.12 / 3.0) 

This query yields 288 objects. After visual inspection we added 32 additional cE candidates to our initial sample. 


S3.3 Estimated effective radius sample 

We notice that well known cE galaxies Messier 32, NGC 4486B, NGC 5846cE have effective radii smaller than 0.4 kpc. We therefore try to 
visually inspect galaxies that satisfy our initial selection criteria but instead of filtering the input catalogue by the velocity dispersion estimate, we 
require galaxies to have the effective radius estimate below 0.5 kpc, similarly to known cEs. In order to do this, we compute mass-to-light ratios 
for galaxies from the input catalogue using the PEGASE.HR stellar population models constructed using the Salpeter IME (36) and then estimate 
effective radii using the published expression (18) assuming that mass follows light. This filter can be formally expressed as 

(corrmag_NUV - corrmag_r - kcorr_NUV + kcorr_r) > 4.0 && 

(corrmag_g - kcorr_g - 25 - 5 * logl0(luminosityDistance(z, 72.0, 0.3, 0.7))) > -18.7 && 

(petror50_r < 2.0 I I petror50_r / 206.265 * luminosityDistance{z, 72.0, 0.3, 0.7) < 0.7) && 

(corrmag_g - corrmag_r - kcorr_g + kcorr_r - gr_fit) > 0.03 && 
age_ssp > 4000.0 && r_eff < 0.5 && z < 0.08 

This query returned 82 objects from which we visually confirmed two. 

To relax the requirement on having a NUV magnitude, we change the virual radius query in a way we did for the UV-dropout sample. The 
corresponding estimated effective radius selection hence becomes: 

(corrmag_g - kcorr_g - 25 - 5 * logl0(luminosityDistance (z, 72.0, 0.3, 0.7))) -18.7 && 

(petror50_r < 2.0 I I petror50_r / 206.265 * luminosityDistance{z, 72.0, 0.3, 0.7) < 0.7) && 

age_ssp > 4000.0 && (corrmag_NUV < -100 I I NOLL_corrmag_NOV) && 

(corrmag_g - corrmag_r - kcorr_g + kcorr_r > 0.8 - 

(corrmag_r - kcorr_r - 25 - 5 * loglO(luminosityDistance(z, 12.Of 0.3, 0.7)) + 20.0) * 0.12 / 3.0) && 

r_eff_virial < 0.5 && z < 0.08 

This query returns 102 objects. Five of them are confirmed visually as cE candidates and included in the preliminary list. 


We attempted to use the effective radii com¬ 
puted by the two-dimensional fitting of SDSS 
images with the GIM2 d software (43). How¬ 
ever, in the vast majority of cases our galaxies 
remain spatially unresolved in SDSS images, 
and the effective radii estimated from the im¬ 
age analysis (43) do not correspond to their 
physical radii. This is illustrated by the fact 
that they are correlated with redshift. There¬ 
fore, we decided not to use these data in our 
further selection and analysis. 

53.4 Additional candidates 

We plot the color-magnitude (g — r, r) dia¬ 
gram and visually inspect all the objects with 
unphysically red colors (also satisfying com¬ 
mon age and velocity dispersion criteria), some 
of which turn out to be galaxies blended with 
bright stars, artefacts and galaxies with under¬ 
estimated g magnitudes by 0.3 mag or more. 
Among the latter group we find 120 cE candi¬ 
dates which we include in the preliminary list 
of cE candidates. 

53.5 Discarding objects with emission lines 

After merging 4 lists of cE candidates selected 
using different approaches described above. 


we obtain one preliminary table with 215 cE 
candidates. We then perform the final qual¬ 
ity check by identifying galaxies with emission 
lines in the spectra from this preliminary list of 
candidates. We compute line fluxes in the spec¬ 
tral fitting residuals and use original flux uncer¬ 
tainties in order to estimate the signal-to-noise 
ratios. Then, we reject the objects where either 
[On] 3730A, or [Oiii] 5007A, or Ha 6565A, 
or [Nil] 6585A exceed the 4-a detection (20 
galaxies). 

Hence, we end up with 195 objects in our 
final sample of cE galaxies. 

S4 Identification of host galaxies 

The reliable formal identification of host galax¬ 
ies is a challenging procedure and is not always 
possible. We employ a manual and to some ex¬ 
tent subjective scenario and note that our host 
identification may be imperfect in cases of ob¬ 
jects with large projected distances from hosts. 

In order to achieve this, we inspect SDSS 
color images in 2 degree vicinities of every 
cE candidate using the CDS ALADIN soft¬ 
ware. In this procedure, we overplot all 
SDSS objects with known spectroscopic red- 
shifts within 2000 km s“^ of the radial veloc¬ 
ity difference from a cE galaxy. Among those 


we visually identify the most luminous galaxy 
in a group/cluster where it is possible. In cases 
where the largest group/cluster member has no 
SDSS redshift (e.g. due to a fiber collision in 
the SDSS spectroscopic targets selection pro¬ 
cedure), we query NASA NED service to ob¬ 
tain its redshift. In a few cases, we select 
a close-by bright galaxy sometimes showing 
signs of interaction with a cE candidate instead 
of the brightest object in the group/cluster. In 
our preliminary list of cE galaxies there are 12 
isolated or “field” objects for which we can¬ 
not identify hosts, whereas in the discarded 
emission line galaxies there is 1 field galaxy. 
Hence, in the final sample we have 11 isolated 
cE galaxies. 

S5 Properties of the final cE sample 

S5.1 Compact elliptical galaxies with tidal streams 

Our final sample includes 8 compact ellipti¬ 
cal galaxies with the visual signs of ongoing 
tidal interaction with their host galaxies (see 
Fig. lS2t resembling the two galaxies with tidal 
streams found earlier (8). The numerical simu¬ 
lations of the tidal stripping (5) demonstrated 
that such features can be observed only dur¬ 
ing a relatively short stage of the cE forma¬ 
tion, that is during the first 200^00 Myr fol- 
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Fig. SI: An example of the full spectral fitting of an SDSS DR7 spectrum of a compact elliptical galaxy with the “NBursts” technique. The black, 
green, red and purple curves represent the observed spectrum, its uncertainties (ilo-), the best-fitting model and residuals. Position of prominent 
atmosphere airglow lines are indicated by dashed blue lines. 


lowing the first close passage of a cE progen¬ 
itor near the host galaxy. Later they quickly 
fade down to very faint surface brightness val¬ 
ues (>29 mag arcsec“^) and become a part of 
the intracluster (intragroup) light. 

S5.2 Host galaxies 

Our new large cE sample allows us to com¬ 
pute statistical properties of the host galaxies 
identified for 184 cEs. Most cEs in our sam¬ 
ple belong to galaxy groups including 5 to 30 
members identified in the SDSS spectroscopic 
sample. Eor 23 cEs, host galaxies are early to 
late-type spirals (in all these cases there is only 
one cE associated to a given host galaxy) while 
for the remaining cases the hosts are massive 
lenticular, elliptical or merging galaxies. 

A distribution of the differences between 
r-band luminosities of cE galaxies and their 
hosts will reflect the distribution of stellar mass 
ratios because stellar r-band M/L ratios of 
cEs and giant early-type galaxies are similar 
(assuming the same stellar initial mass func¬ 
tion). The obtained Mr{cE) — Mr {host) dis¬ 
tribution is shown in Eig. IS3I It is skewed and 
peaks at 3.25 mag corresponding to the stel¬ 
lar mass ratio of 20. Giant elliptical galaxies 
contain about 90% of dark matter within their 
virial radii (44, 45), while for nearby cEs like 
Messier 32 the measured dark matter content is 
negligible being consistent with the tidal strip¬ 
ping scenario where the extended dark matter 
halo is stripped before the stellar component. 
This yields in the median cE-to-host mass ratio 
of about 200. 

Worth mentioning that prototypical cEs, 


Messier 32 and NGC 4486B with the cE-to- 
host stellar mass ratios of 100 and 60 respec¬ 
tively are sitting in the tail of the distribution 
including a few per cent of objects. 

In Eig. lS4l we show the distribution of pro¬ 
jected distances on the sky between cEs and 
their identified host galaxies. The distribu¬ 
tion looks almost uniform between 20 kpc and 
1 Mpc (in logdproj) with a prominent peak 
at 50-100 kpc. The numerical simulations of 
tidal stripping (5) demonstrate that if a cE pro¬ 
genitor comes too close to the host galaxy, the 
quick orbital decay causes its very quick ac¬ 
cretion on-to the host on a timescale of a few 
hundred Myr. On the other hand, under typi¬ 
cal conditions in clusters and rich groups, tidal 
stripping becomes significantly less efficient at 
distances exceeding 200 kpc, and this explains 
the decline. The high-dproj tail probably con¬ 
tains mostly cEs which were gravitationally 
ejected from the cluster/group centers by three- 
body interactions. 

S6 Construction of caustics for cE host clusters 
and groups 

We identify 191 cE galaxies as members of 
galaxy clusters from the groups and clusters 
catalogue of (20). 7 of them have the luminos¬ 
ity rank 1 in a group of 2 galaxies and, there¬ 
fore, were considered isolated or field cEs. In 
order to study the dynamical connection be¬ 
tween cE galaxies and their host environments, 
we plot phase space diagrams showing red- 
shifts of clusters members versus their pro¬ 
jected distances from cluster centers. Such 
diagrams usually show a characteristic trum¬ 


pet shape while their boundaries are called 
caustics. (23) provided a technique to compute 
them under assumptions of spherical symmetry 
and hierarchical clustering for the formation of 
the large-scale structure. This method works 
as a convenient way to estimate the escape ve¬ 
locity from the gravitational potential well gen¬ 
erated by the cluster at each cluster-centric ra¬ 
dius. We used the implementation of (23) tech¬ 
nique in the Caustic App version 1.2 (46) to 
find caustics in the data for our galaxy clus¬ 
ters and identify cE positions with respect to 
the caustics. 

We also create an ensemble galaxy cluster 
to compare the distribution of normal cluster 
members to the one of our cE candidates in 
the phase space. We take 33 galaxy clusters 
that have more than 20 members and include 
at least one cE candidate. Their masses, ve¬ 
locity dispersions and effective radii were al¬ 
ready known after the previous step when their 
individual caustics were computed and anal¬ 
ysed. We project every cluster on a tangen¬ 
tial plane with the projection center at the clus¬ 
ter center, and then normalize projected coor¬ 
dinates of galaxies by half-mass radii of the 
corresponding clusters. The relative velocities 
with respect to a cluster center are normalized 
by the cluster velocity dispersions. Then, all 
the clusters are shifted to the redshift Zens ~ 
0.05. In this fashion, positions and radial ve¬ 
locities of 2,592 galaxies in 33 galaxy clusters 
are brought to the same scale ready to be used 
as ensemble cluster. We computed caustics of 
generated ensemble cluster using default set¬ 
tings of Caustic App. 
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Fig. S2: Compact elliptical galaxies showing signs of ongoing tidal interactions with their hosts. 

ce.objid 

ra 

dec plate 

mjd fiberid 



587722953304375642 237.29089 0.21562 342 51691 548 


Tab. SI: Properties of 195 cEs and their host galaxies found in this study. Missing value in the host.objid column indicates that this row 
corresponds to an isolated cE galaxy without an identified host. The complete table along with column descriptions is available in Microsoft Excel 
format as a separate supplementary file. 



Fig. S3: A distribution of r-band magnitude difference of cEs and 
their identified host galaxies M^e - Mhost reflecting the stellar 
mass ratios. 
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